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We focus on the features of spectra and diagrams of states obtained via exact diagonalization technique for 
finite ionic conductor chain in periodic boundary conditions. One dimensional ionic conductor is described 
with the lattice model where ions are treated within the framework of "mixed" Pauli statistics. The ion transfer 
and nearest-neighbour interaction between ions are taken into account. The spectral densities and diagrams 
of states for various temperatures and values of interaction are obtained. The conditions of transition from 
uniform (Mott insulator) to the modulated (charge density wave state) through the superfluid-like state (similar 
to the state with the Bose-Einstein condensation observed in hard-core boson models) are analyzed. 
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1. Introduction 

Ionic conductors are a wide class of physical and biological objects ranging from ice to DNA 
membranes. One of the most interesting subclasses of these are superionic conductors that exhibit 
high temperature phase with high conductivity that arises due to the motion of ions [lj or pro- 
tons 0. Theoretical description of systems with ionic conductivity is most frequently based on 
the lattice models. Some of them treat ions as Fermi-particles focusing on different aspects of the 
ionic subsystem like long-range interactions or interaction with phonons 0, 01. Some recent 
attempts have also been made towards short-range interactions between particles [8l4l2 1 . 

However, a more correct consideration of ions should be based on the mixed statistics of 
Pauli [l3[ since these particles are bosons by nature but they also obey the Fermi rule. Due 
to the special commutation rules, the utilization of Pauli operators generates additional mathe- 
matical complexities. On the other hand, this approach might be very effective. For instance, it 
has been shown that the lattice model of Pauli particles is capable of describing the appearance 
of superfluid-like state (that corres pon ds to superionic phase) in the system even in the absence 
of interaction between particles |l4l4l6j . On the other hand, the lattice model of Pauli particles is 
similar to the hardcore Bose-Hubbard model widely used for the description of ionic conductivity 
phenomena as well as for the modelling of energy spectrum of absorbed ions on a crystal surface 



and intercalation in crystals [17(. Bose-Hubbard model also exhibits the transition from Mott in 



sulator state to superfluid-like state H|-24|. Some of the authors also observe the possibility of 
formation of intermediate "supersolid" phase that may appear on the phase diagrams alongside the 
transition from dielectric (CDW) to superfiuid phase. 

In this work we focus on the diagrams of state for one-dimensional ionic conductor described 
by the system of Pauli particles. Our lattice model includes ion transfer as well as the interaction 
between nearest-neighbouring ions. We calculate the single-particle spectral densities of the finite 
system in periodic boundary conditions and obtain the diagrams of state analyzing the features 
of this spectra. The conditions of transition from Mott insulator (Ml)-like state to the modulated 
charge density wave (CDW) state through the superfluid(SF)-like state (similar to the state with 
the Bose-Einstein (BE) condensation observed in hard-core boson models) are discussed. 
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2. The model for ionic conductor 

Let us consider the chain of heavy immobile ionic groups (large circles in figure 1) and light ions 
that move along this chain occupying positions denoted by small circles in figure 1. The subsystem 
of light ions is described with the following Hamiltonian 

H = ty](cfcj + i + cf +1 Cj) + Vy^njm+i - fiy]nj. (2.1) 

i i i 

This model takes into account the nearest-neighbour ion transfer (with hopping parameter t) and 
interaction between ions that occupy nearest-neighbouring positions (with corresponding parame- 
ter V) . If this Hamiltonian is considered within the framework of Fermi statistics, the corresponding 
model is known as spinless-fermion model. This model is widely used in the theory of strongly cor- 
related electron systems [25j as well as for the description of ionic conductors [2|| . A more complex 
two-sublattice case of this model can be applied to a proton conductor [27j ■ More correct consid- 
eration of ions should be based on "mixed" Pauli statistics and this approach is used onwards. In 
this case, the model (|2.1|) is equivalent to the extended hard-core boson model, i.e., boson Hubbard 
model with repulsive interaction between nearest neighbours and infinite on-site repulsion [28]. The 
latter is often applied to the investigation of the problems of BE-condensation and superfluidity. 




Figure 1 . The model for one-dimensional ionic conductor. Large circles denote heavy ionic groups 
while the small ones denote light movable ions. 



3. Exact diagonalization technique 

We calculate the spectral densities of one-dimensional ionic Pauli conductor using exact diag- 
onalization technique. For the chain of N sites, we introduce the many-particle states 

| ni, a n^b . . . n N . a n Nib ). (3.1) 

The Hamiltonian matrix on the basis of these states is the matrix of the order 2 N x 2 N and is 
constructed as follows: 

N 

H mn =J2[t (H£l + H%1) + VH& - , (3.2) 

»=i 

H rnl = (ni...\c+c l+1 \n' 1 ...) = 6{n l -n' t -l)S(n l+1 -n' l+1 + l) 

x n 5 ( n ' _n 0) 

H ml = (ni...\cf +1 Ci\n[...) = S(n i -n' i +l)S(ni +1 -n' i+1 -l) 

n 5 ( n i- n i)i 



where 



X 

H ml = (ni...\nin i+ i\n[...) = S(ni-l)6(n' i -l)d(n i+ x-l) 
x<5«+i-l) I] *(»i-nD. 

HW = (n 1 ...\n t \n' 1 ...)=6(n l -l)6(n' l -l)l[6(n l -n' l ). 



23702-2 



Spectral densities and diagrams of states of 1 D ionic Pauli conductor 



This matrix is diagonalized numerically 

U~ l HU = H = ^2\ p X pp , 



(3.3) 



where X p are eigenvalues of the Hamiltonian, X pp are Hubbard-operators. The same transformation 
is applied to the creation and annihilation operators 



U- X CiV = J2 A V1 Xm , U^C+U = A rs XrS 



(3.4) 



which are required to construct one-particle Green's function <C Ci, a \cj a 3> that contains infor- 
mation about one-particle energy spectrum of the system. For Pauli creation and annihilation 
operators, this Green's function can be constructed in two ways, i.e.. commutator Green's function 



« Cl (t)\ci(t') »= -ie(t-f)<[ci(t) 1 q + (i')]> 

and anticommutator Green's function 

« a(t)\4(f) »= -«*(* - 0({ci(*),c?-(0}>- 



(3.5) 



(3.6) 



Imaginary part of these Green's functions are one-particle spectral densities (also referred to as 
densities of states or DOS) 



1 N 

= --»f X] 1x11 <C c ^\ c ta > 



ttN 
1 



i=l 

N 



V A A* ' 



rje' 



U-(Xq- Xp ) 



(3.7) 



where Z = Ylp e ^^ p - Spectral densities in (|3.7|l . obtained from commutator rj = 1 (|3.5|) and 
anticommutator r/ = —1 (|3.6p Green's functions, respectively, exhibit a discrete structure, i.e., 
consist of several <5-peaks due to the finite size of a cluster. Therefore, we apply the periodic 
boundary conditions to the cluster and introduce a small parameter A to broaden the ^-peaks 
according to Lorentz distribution 

^b^- (3 - 8) 



4. Results and discussion 

We perform calculations of one-particle spectral densities of one-dimensional ionic Pauli con- 
ductor (|2.1| for the chain of ten sites (N = 10) in periodic boundary conditions. To test the results 
we compare them with the exact solution obtained by means of fermionization procedure fl4l | in 
the absence of nearest neighbour interaction (V = 0) for different values of chemical potential (fig- 
ure [2]) . The redistribution of statistical weight with the change of chemical potential level can be 
observed, and a good level of agreement is achieved. The analysis of spectral densities is a way to 
distinguish the different states of an ionic subsystem and the corresponding conditions. In the case 
of half-filling, as we turn on the interaction starting with spectral density whose shape corresponds 
to SF state, we observe the development of the gap in the spectra (figure [3]). A similar effect was 
found for ionic and proton conductors described by the similar models within the framework of 
Fermi statistics. It was connected with the splitting of spectra due to charge ordering with doubling 
of the lattice period (ill [l2| . The detailed analysis of spectral densities of Fermi and Pauli models 
of ionic conductor can be found in j29j. So, here we have a transition from SF to CDW state. This 
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Figure 2. Commutator spectral densities of non-interacting (V = 0) ionic Pauli conductor for 
different values of chemical potential (right figure, e = — (x) compared to exact results obtained 
in [14 | via fermionization procedure (left figure, e = eo — fJ-)- t = 1,T = 0.2, A = 0.4. Spectral 
density on the left figure is scaled to 2tv, while on the right figure it is scaled to unity. 




Figure 3. Commutator spectral density of ionic Pauli conductor for different values of nearest- 
neighbour interaction V at half filling (n = 0). t = 1, T = 0.2, A = 0.4. 



transition was described analytically for T — in [3(1 131 1 . We have investigated the case of T ^ 
when this transition manifests itself as a crossover. 

The existence of a gap in the spectra which separates the bottom of the energy band from the 
chemical potential level is also the sign of the presence of homogenous MI state. The vanishing 
of the gap and the appearance of a negative branch points to the transition to SF-like state 
(see, for example, [32]). According to these criteria we analyze the spectral densities at different 
temperatures and values of interaction and build the corresponding diagrams of state (figure 0]) . 
The system is in MI homogenous state at high temperatures and far away from half-filling (at 
large $). As the temperature decreases or one comes closer to half- filled case, the system undergoes 
transition to SF-like state which corresponds to the appearance of the negative branch without any 
gap on the spectral density. At further temperature decrease and closer to half filling, we observe a 
transition to the state with the gap on the spectral density that corresponds to the CDW-ordering 
(the negative branch still exists). As the interaction strength V increases, such a region becomes 
broader while the region of SF-like state becomes smaller. On the other hand, with the decrease 
of V, the CDW state diminishes and disappears at V « t. It should be mentioned that the sequence 
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of states that the system is going through at the increase of mean occupancy 6, corresponds to the 
phase diagram obtained in [3CJ, [3jJ, |33j . We have also performed a detailed analysis of the transition 
to SF-like state at different temperatures and values of interaction (figure [5]) . It is interesting that 
at weak interactions (V < 1), the increase of interaction strength facilitates the formation of 
SF-like state while further increase of V suppresses this transition. The shift of the curves that 
separate MI- and SF-states towards smaller values of 8 with the increase of interaction strength 



that we have obtained is also observed on the phase diagrams obtained by other authors |30L 131 1 
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Figure 4. Diagrams of state for different interactions and the spectral densities for negative 8 
that correspond to CDW, SF and MI states, t = 1, A = 0.25. 8 = (n) — 1/2 denotes the deviation 
from half-filling. 
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Figure 5. The curves that separate SF- and MI- states at different temperatures and values of 
interaction. SF-state is on the left hand side of the curves and below them, t = 1, A = 0.25. 
S = (n) — 1/2 denotes the deviation from half-filling. 



5. Conclusions 

We have performed the analysis of diagrams of states of one-dimensional Pauli ionic conductor 
using an exact diagonalization technique. We have shown that the system undergoes transition 
from Mott insulator to superfiuid-like state and then to CDW-sate. At weak interaction, the latter 
transition may vanish. 
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CneKTpa/ibHi rycTi/im/i Ta fliarpaivin CTaHy OflHOBi/iMipi-ioro ioHHoro 
npoBiflHi/iKa l~lay/ii 

I.B. CTactoK, 0. Bopo6£ioB, P.9\. CTei_n'B 

I HCTMTyT 4>i3HKH KOHfleHCOBa h wx cwcTeM HAH y Kpa ihi/i , By.n. I. Cb6h L4i L4bKoro, 1, 79011 JlbBiB, YKpaiHa 

Po6oTa npucBaneHa bh BHeHHio eHepreTi/iHHoro cneKTpy Ta fliarpaM CTa HiB, OTpwMaHnx MeTOflOM tohhoT flia- 
roHa/ii3aL4iT fl/ia CKi h neHHoro ioHHoro .na h L4K)roBoro npoBiflm/iKa b nepioflWHHWx rpa h n h h hx yMOBax. OflHOBi/i- 
MipHM m ioHHM m npoBiflHUK oni/icyeTbcn rpaTKOBOK) MOfleji/iio, fle ioHW po3r/i;RflaK)TbCiR aK nacTUHKH riay/ii, npi/i 
LibOMy BpaxoByeTbca ioHHufi nepeHoc i flBonacTHHKOBa B3aeMOfliiR Mi>K Haw6.nn>KHHMH cyciflaMW. Ey.no po3pa- 
xoBaHO Ta npoaHajii30BaHO cneKTpajibHi rycTHHH Ta fliarpaMW CTaHy TaKoT cucTeMW fljia pi3Hwx TeMnepaTyp 
Ta Be/i m h m h B3aeMOfliT. ripoaHa.ni30BaHO y mobh nepexofly ci/icTeMU 3 oflHopiflHoro (cTaHy t. 3B. MOTTiBCbKoro 
flie/ieKTpuKa) y MOfly/ibOBaHMM CTa h Hepe3 CTaH Ti/iny cpa3H 3 6o3e-KOHfleHcaTOM (nofli6HoT flo Hafln/inHHoT 

Cpa3W B MOfle/iax >KOpCTKWX 6030HiB). 

KyiKJHOBi enema: CTaTtACTWKa llayjii, rycTMHa CTamB, iohhiaia npoBiRHMK 
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